Correction to the Casimir force due to the anomalous skin effect 
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Abstract 



The surface impedance approach is discussed in connection with the precise calculation of the 

CO 

Casimir force between metallic plates. It allows to take into account the nonlocal connection 
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^ . be described by two impedances Z s (u>,q) and Z p (uj,q) corresponding to two different polarization 



between the current density and electric field inside of metals. In general, a material has to 
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the Thomas-Fermi screening. The nonlocal correction to the Casimir force is calculated at zero 
temperature. This correction is small but observable at small separations between bodies. The 
same theory can be used to find more significant nonlocal contribution at uj ~ ui p due to the 
plasmon excitation. 



states. In contrast with the approximate Leontovich impedance they depend not only on frequency 
uj but also on the wave vector along the plate q. In this paper only the nonlocal effects happening at 
frequencies uj < u p (plasma frequency) are analyzed. We refer to all of them as the anomalous skin 
effect. The impedances are calculated for the propagating and evanescent fields in the Boltzmann 
approximation. It is found that Z p significantly deviates from the local impedance as a result of 
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I. INTRODUCTION 

The Casimir force |j] between uncharged metallic plates, predicted in 1948 as a quantum 
electrodynamics effect, only recently became a subject of systematic experimental investi- 
gation. The reason is that nowadays with the development of microtechnologies a reliable 
control of the separations between bodies smaller than 1 [im became possible. A variety of 
methods have been used to measure the force. In the torsion pendulum experiment |2|, first 
in the modern series, the force between a sphere and a plate covered with gold was measured 
with the accuracy of 5%. A significant progress was achieved in the atomic force microscope 
(AFM) experiments 0, Q], where the sphere was attached to a cantilever. The force was 
deduced from the cantilever bending when the plate was approaching the sphere. In this 



4, 



experiment the force was measured with a 1% precision. The best result using AFM 
was found when the sphere and the plate were covered with gold and special care was taken 
to control the surface roughness of metal. The same precision was reached in the crossed 
cylinders experiment jsj], where extremely smooth gold films were used. Sophisticated micro- 
electromechanical system (MEMS) [6] was used to measure the force between the gold plated 
sphere and a suspended paddle. It demonstrated the nonlinear behavior of the mechanical 
oscillator due to the Casimir force. The only experiment, where the force was measured in 
the plate-plate geometry [7(, was done using an oscillating beam whose resonance frequency 
changed in response to the force. Up to date, the most precise experiment j^, |t| explored 
the MEMS device similar to that in Ref. p]. The precision was improved due to the use 
of the dynamical method. Additionally the change in the resonance frequency of the me- 
chanical oscillator was measured using the phase jump instead of resonance behavior of the 
amplitude. In this way the force was found with the relative accuracy of 0.25% |sl IqL llo|. 

To draw any conclusion from the experiments one has to predict the force theoretically 
with the precision comparable with the experimental errors. In its original form, the Casimir 
force [ill 



= (i) 



7T 2 He 

240^ 

was calculated between the ideal metals. It depends only on the fundamental constants 
and the distance between the plates a. The force between real materials, described by its 
dielectric functions e (u), was deduced for th. flrst to by L ifsUtz Q. Co.eeUoee t o 
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Eq. (P) can be quite significant at small separation between bodies. To calculate the force 
with hig h precision the Lifshitz formula is used with the optical data taken from handbooks 
[13I [lj The data are available only up to some low cut-off frequency u cut . For good 
metals such as Au, Al, Cu the data can be extrapolated to lower frequencies with the Drude 
dielectric function 



e(w) = l 



LOt 



(2) 



uj (oj + iuj T ) ' 

which includes two parameters: the plasma frequency u p and the relaxation frequency uj t . 
These parameters can be extracted from the optical data at the lowest accessible frequencies. 
In this way the force has been calculated [la, 0, 

with the highest possible 

precision. There is some disagreement between the results of different authors connected 



with the choice of the relaxation frequency ui T jis| . 



This frequency can be found by fitting the 



ow frequency optical data with Eq. (J2J) [l7[ or extracted from the bulk material resistivity 
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. These detai 
In Refs. Q, 



s are important at small separations. 

it was stressed that using the handbook optical data one finds not 
the actual force but rather the upper limit on the Casimir force. The reason is that the 
handbooks comprise the data for the best samples; any material imperfection will reduce 
the reflectivity and, as a result, the force will be smaller. In the experiments the force is 
measured between the bodies covered with a metal. The metal is deposited on a substrate 
with the evaporation or sputtering technics The resulting film thickness is typically 

in the range 100 — 200 nm. It was already noted 0, 13 that the optical properties of the 
films can deviate significantly from those of the bulk material. The main reasons for these 
deviations such as voids in the films and electron scattering on the grain boundaries were 



indicated recent 
and MEMS 
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2 it ] and the methods to estimate the effects were outlined. For the AFM 



experiments the influence of these effects on the force was estimated on 



the level of 2%. Additional work from the theoretical and experimental sides has to be done 
to refine these estimates. 

An alternative way to calculate the Casimir force using the surface impedance of metals 
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24] in the literature. The general 



instead of the dielectric function was discussed 
formula for the Casimir force in this approach was deduced for the first time in Ref . [2a] . It is 
the same Lifshitz formula but the reflection coefficients are expressed in terms of the surface 

n 

impedance. In the cited papers the approximate Leontovich impedance 26] depending only 



on frequency was used. The hope was that in this approach it will be possib 



the lon g st anding problem with the temperature correction to the force 
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32J, |33j, |34|. However, although the impedance approach itself seems reasonable and 
well motivated, the use of the approximate impedance for the Casimir force happened to be 
unjustified. The Leontovich impedance is well suited for the propagating photons impinging 
on a metal but for the Casimir problem the exact impedance has to be used since the 
important contribution in the force comes from the evanescent electromagnetic field. The 
exact impedance was shown to give the same result for the force as that in the dielectric 
function approach [3^ . 

The surface impedance is the only way to describe the interaction between the electro- 
magnetic field and metal in the case when the relation between current and electric field 
in the metal becomes nonlocal. For example, at low temperatures the mean free path of 
the electrons in metals becomes larger than the field penetration dep th, and the relation 
between the current and field becomes nonlocal at low frequencies |3Jj, |3j| and the anoma- 
lous skin effect is realized. In connection with the Casimir force the contribution of this 
effect in the temperature correction was discussed in Ref. [23 1. On the other hand at high 
frequencies u ~ u p the charge density fluctuations can propagate in the material (plasmons) 
taking away the energy from the incoming field. This is also an example of nonlocal effect 

38j | at small separa- 



which was shown to give a significant correction to the Casimir force 
tions between bodies. In the nonlocal case both time and space dispersion happen and the 
dielectric function depends not only on frequency u but also on the wave vector k. Actually 
in this case one has to separate two dielectric functions: the longitudinal function ei (u, k), 
which describes material response to the longitudinal (in respect to k) electric field, and 
the transverse dielectric function e% (u, k) describing the response to the transverse field. A 
systematic way to calculate the surface impedances via the dielectric functions and explicit 
expressions for these functions were given in a series of classical papers by Kliewer and Fuchs 

0yQ. 

The aim of this paper is to provide the basis for systematic investigation of the corrections 
to the Casimir force due to the nonlocal effects. Here we consider only the minor corrections 
which appear at low frequencies uo « u p but introduce a general approach which is true at 
any frequencies. This approach is not new in the condensed matter physics but it has never 
been discussed before in connection with the Casimir effect. 
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We deliberately do not consider here the correction to the force in the nonzero temper- 
ature case though our conclusions about the low frequency behavior of the impedances for 
two polarization states will be important for the discussion of the temperature correction. 

The paper is organized as follows. In Sec. Othe main definitions are introduced allowing 
to calculate the impedances for two polarization states via the longitudinal and transverse 
dielectric functions. The explicit expressions for the nonlocal dielectric functions in the 
Boltzmann approximation are given. The impedances are calculated first for propagating 
fields to compare the calculations with the known results. Then we discuss the continua- 
tion procedure to the range of evanescent fields and calculate the impedances at imaginary 
frequencies to put them later in the Lifshitz formula. The low frequency behavior of these 
impedances important for the thermal correction to the Casimir force is discussed specif- 
ically. In Sec. IIHI the actual calculations of the correction to the Casimir force due to 
nonlocal effects are presented for plate-plate and sphere-plate geometry. The discussion and 
concluding remarks are given in Sees. IIVI and Ivl 



II. SURFACE IMPEDANCES OF METALS 



The anomalous skin effect was incorporated into the general theory of the optical prop- 
erties of metals with the detailed paper by Reuter and Sondheimer j^j]; the qualitative 
description was given earlier by Pippard 43]. Pippard was the first to point out that, in 
general, the electric field inside a metal cannot be considered as spatially constant. In gen- 
eral, the connection between the current and the field becomes nonlocal. The current is 
given by a definite integral involving the values of the electric field at all points in the metal, 
and the Maxwell's equations therefore lead to integro-differential equation from which the 
electric field has to be determined. The expression for the current was deduced at the 
conditions that the electrons can be considered as quasi-free and the collision mechanism 
can always be described in terms of a mean free path I or, alternatively, a relaxation fre- 
quency u T . The mean free path was assumed to be independent on the direction of motion. 
An additional assumption about the electron reflection off the surface was introduced by 
Pippard: a fraction P of the electrons arriving at the surface is scattered specularly, while 
the rest are scattered diffusely. Only normal incidence of the electromagnetic field on the 
metal was considered in Ref. Q], however, it was stressed in the Kliewer and Fuchs paper 
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3^1 that the theory of anomalous skin effect can be considered complete only when arbitrary 
incidence is taken into account. 

The relation between the current density and the electric field can be found solving the 
Boltzmann equation if the free electrons and field wavelengths large compared to the wave- 
length of an electron at the Fermi surface. For the description of the anomalous skin effect it 



is a valid approximation. In genera 
using the linear response theory [44 1 



, this relation can be found in different approximations 



A. Propagating waves 

Following Kliever and Fuchs [39] we consider a plane wave of angular frequency uo in- 
cident from vacuum at an angle •& upon the surface of metallic half-space. The geometry, 
together with the choice of the coordinate system, is shown in Fig. ^ One can separate 
two polarization states for the wave. For the s polarized wave the electric field is directed 
in the y axis, while for the p polarized wave the electric field is in the x — z plane and 
has nonzero z component. For clarity let us sketch out how the specific expressions for the 
surface impedances have been deduced in Ref. j^. 

For the s polarized wave the field can be written in the form E = E y (z) exp(— iut + 
ik x x)n y , where is the unit vector along the y axis and k x = (uj/c) sin $ is the ^-component 
of the wave vector in the incoming wave. Since this component will play significant role in 
what follows, we will use for it also a special notation q = k x , which is settled in the field of 
the Casimir force. The z dependence of the electric field E = (0, E y , 0) for the s wave inside 
of metal can be described with the Maxwell equation 

^ ~ klE y + ^D y = 0, (3) 

where D is the displacement field. This equation is valid for z > 0. Taking the Fourier 
transform of Eq. (jSJ) over the z coordinate, we obtain 

- « + kl)e, + p y = -MH!) _ ^L2>, (4) 

where S and T> are the Fourier transformed fields defined as 
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S(k z ) = J dzE(z)e- tkzZ (5) 

— oo 

and similarly for T>. Equation (J3J) describes the behavior of the electric field in an infinitely 
extended medium. Furthermore, the right hand side in Eq. (@J) is undefined until we find 
a relation between the two derivatives at the surface. This relation involves describing or 
modelling the surface. In this work we assume that the electrons scatter elastically from 
the surface (specular reflection). This assumption is equivalent to assuming an infinitely 
extended medium, needed to obtain Eq. (0J), since an electron bouncing from the surface 
cannot be distinguished from an electron coming from a fictitious medium on the vacuum 
side. This is taken into account imposing the symmetry requirements 



E y (z) = E y (-z), D y (z) = D y (-z). 

Thus, using the Maxwell equation 

dE y (+0) dE y (-0) 



dz 



dz 



-i (u/ c) H x 



(6) 



(7) 



and from Eq. (0J) one finds 



2iuj 



k 



k -\- k 



(8) 



H x (+0) c (LU 2 /c 2 )s t -k 2 

The inverse Fourier transform of this equation evaluated at z = gives the surface impedance 
for s polarization: 



Z s (k x , u) = 



_ E y (+0) _ i u f dk z 

7T c J (lo 2 /c 2 ) e t — k 2 



(9) 



— oo 



H x (+0) 

For p polarization the problem is slightly more complicated since we have two non- vanishing 
components of the electric field E x and E z . Following a similar line of reasonings as before, 
the impedance for p polarization is obtained as: 



Z p (k x , u) 



E x (+0) 



oo 

i u I dk 7 



+ 



(10) 



H y (+0) 7T c J k 2 [(uj 2 /c 2 )ei (io 2 /c 2 )e t -k 2 _ ' 

— oo 

It is natural that both the longitudinal e\ (k, iS) and transverse e t (k, uo) dielectric functions 
contribute to Z p because in the p wave the electric field has both components. 



Since the imp edance approach caused recently some confusion in the field of the Casimir 
22], |23j, |2^ , a few comments concerning the impedances Eq. © and Eq. (jlUj) should be 



force 

made. First, there is not one but two impedances corresponding two different polarizations. 
Second, the impedances depend not only on the frequency but also on the wave vector along 
the metal surface q = k x = (id / c) sin$. Only for the normal incidence the impedances for p 
and s polarized waves coincide and depend only on frequency. Third, no specific assumptions 
about the dielectric functions E\ and et were made in the derivation of Eq. Q and Eq. (JTUJ). 
In particular, the local functions e% (0, u) = e t (0, id) — e (id) can be used. In this case the 
integrals can be easily calculated to find so called classical or local impedances 
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\je{u)- (cq/u)' 



E(ld 



y/e{u)- {cq/u) 2 

These expressions completely coincide with those introduced in Refs. j^, 



(11) 



38j and, as was 



shown there, exactly reproduce the Casimir 



approximate one 



22 



23, 



brce in the dielectric function approach. The 



24j from the beginning was introduced as the 



Leontovich impedance used in Refs. 

nn 

36l . l37| for applications in radiophysics. For the propagating waves it 



really has sense because cq/u < 1 but for metals in the microwave range \e (u)\ ^> 1. So 
one can neglect the wave vector along the plates (q = k x ) in the impedances to get just one 
frequency dependent function. However, in the Casimir force significant contribution comes 
from the evanescent fields for which cq/u > 1. In this case the Leotovich approximation 
fails especially in the limit uo — > 0, which is important for the analysis of the temperature 
correction. 

The dielectric functions were found (3^ solving the Boltzmann equation and the result 
is the following 



e t (k,uj) = l + X iB H - 



-ft («) 



(12) 



ei (k,u) = 1 + Xib H - 



uo: 



:fi («) 



(13) 



UJ (id + lid T ) 

where the phenomenological susceptibility xib (^) was introduced to describe the interband 
transitions since these processes are beyond the quasi-free electron model. The functions 
ftj (u) taking into account nonlocal effects are defined in the following way 



S 



/«(«) 



3 

2^ 



u 



(14) 



1 n + w 

-ii H — m 

2 Vl-u 



.u; r % oj T ^ 



where the variable u responsible for the nonlocal effects is 



1 + u 



u 



(15) 



u = — — — , 16 

UJ + 10J T 

and Vp is the Fermi velocity. In the local limit k — > both functions (|T2*|) and (fT3*j) reduce 
to the local dielectric function (Drude plus interband transitions) 



e (u) = 1 + xib (w) - 



(17) 



The function e t was found first by Reuter and Sondheimer |42j]. Since these authors 
considered only the normal incidence, in their work there was no This function appears 
at non- normal incidence. In the p wave there is the normal field component E z giving rise 
to charge fluctuations to which the system responds via the longitudinal dielectric function. 
It should be mentioned that with the charge fluctuations the relaxation of the perturbed 
electron distribution toward equilibrium was chosen to the local state of charge imbalance 
but not to the uniform distribution. The denominator in Eq. (|T3j) describes this effect. 

The equations (fI2 |) -([T6' J) are used in the optics of metals 45( to predict the reflectance or 
absorptance of the materials. It is easy to find the reflection amplitudes r s and r p for s and 
p polarizations expressed via the impedances as: 



u> 2 
7T 



Q 



The reflectance and absorptance are given by: 



CJ 2 
7T 



-q 2 + - C Z P 



(18) 



R 



s,p 



A, 



s,p\ 



(19) 



In what follows we will use dimensionless variables, which are more convenient for nu- 



merical calculations. We define 
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To verify the procedure we recalculated the absorptance with 7 = and the Fermi 

velocity Vp = 0.85 • 10 8 cm/s (potassium) to compare with the same calculations in Ref. 

. The results are presented in Figs. El and El Absorptance at the normal incidence d = 
(Q = 0) is shown in FigEl In this case both polarizations give the same result. The nonlocal 
case is presented by the solid line. The absorptance in the local limit calculated with the 
impedances is given by the dashed line. The usual increase in the absorptance can be 
seen at low frequencies Q ~ 10~ 3 due to the anomalous skin effect. For the incidence angle 
$ = 75° the absorptance of the p polarized wave is shown in Fig. El In this case there 
is an additional peak in absorptance at higher frequencies SI ~ 0.1. It appears only for p 
polarization; the s polarization behaves similar to the case = 0. At smaller 7 both of 
the peaks become much more significant. These results are in full agreement with those of 



Kliewer and Fuchs 



B. Evanescent fields 

The fluctuating currents in the plates are the sources of fluctuating electromagnetic fields 
responsible for the Casimir force. The typical separation between bodies in the Casimir 
force experiments is smaller than the wavelength A of visible light. If we consider one plate 
as an emitter and the other one as a receiver, then for significant part of the spectrum 
contributing in the force the receiver will be in the near field zone of the emitter. In this 
case the propagating field radiated by the emitter will be small in comparison with the 
evanescent field which exists around the emitter at the distances ~ A. The well known 
example of such an emitter is the Hertz dipole. At small distances from the dipole ujr/c <C 1 
one can neglect the retardation and the field around the source is just the field of the static 
dipole decaying as 1/r 2 . When the force is calculated using the Green function method jl^ . 
the Green function is exactly the dipole field modified by the presence of the plates. The 
planar geometry of the problem makes it preferable to expand the dipole field on the plane 
waves. The plane waves obeying the condition u 2 /c 2 < q 2 do not propagate in the gap 
because the normal component of the wave vector is pure imaginary. 



There were some speculations in the literature (see, for example, 2M )> inspired by the 
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problem with the temperature correction to the Casimir force, that for the evanescent fields, 
the standard expressions for the Fresnel reflection coefficients should be modified. In this 
connection we have to stress that the evanescent fields are the subject of the near-field 
optics 46| (see also [47[ for a review), where standard electrodynamic approaches are used. 
Additionally, the longitudinal dielectric function can be probed in the evanescent range by 



the scattering of a beam of charged particles or fast electrons from the material 48] . In this 
way the function Im{l/ei(k, uj)} can be directly extracted from the experiment, where k is 
connected with the momentum and u with the energy losses of the charged particles. No 
necessity for modification of the standard electrodynamics was noted so far. A consistent 
way for the description of evanescent fields is just the analytic continuation of the Eqs. Q, 
(ITU1) . (H2J)-(IinD, and Eq. (HEJ) to the range u 2 /c 2 < q 2 . 

Originally the Lifshitz formula for the Casimir force was written as an integral over real 
frequencies u h|. In this representation one has to calculate first the integral over the 



variable p = a/1 — (cq/oo) 2 in the range < p < 1 (propagating fields) and then integrate 
over the imaginary axis p — i \p\ from zero to infinity (evanescent fields). So the propagating 
and evanescent fields were clearly separated. The alternative representation of the same 
formula [12] is more popular because of faster convergence of the integrals. In this case the 
integration is done over the imaginary frequencies uj = i( but the inner integral over p = 



\/l + (cq/() 2 is calculated from 1 to oo. Formally we are always in the evanescent domain 

because at imaginary frequencies the normal component of the wave vector is pure imaginary 

k z = i^/C 2 /c 2 + q 2 . For this reason we will not investigate specially the domain q 2 > uj 2 /c 2 

making the analytic continuation on q but instead we will make the analytic continuation 

to imaginary frequencies. This procedure is well defined for the response functions which 

are analytical in the upper half of complex plane uj. In the electrodynamics the response 

2i-i 



functions are the components of the Green 
but not the dielectric functions themselves 



unction e l 1 (k, uj) and [(uj 2 /c 2 ) e t (k, uj) — k 



Exactly these expressions take part in 



the impedances Q and (fTU|) and, therefore, the impedances can be considered as analytical 
functions of uj. 

Using the dimensional variables (J2U|) the impedances at imaginary frequencies (Q — > ifl) 
can be written as 
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2Vt C cosh x dx 
kQJ cosh 2 x + §k 



z s (Q,n) = ~l l2 g 2 x ^ x , (2i) 





oo 



z P (Q,n) 2<J f ' lx 



7r Q J cosh x 







+ 



cosh 2 x — 1 



(22) 



(fl, v) cosh 2 x + §kt (fi, ^ 
Here we introduced a new variable of integration x which is defined by the relation k z = 
k x sinh x- For the dielectric functions at imaginary frequencies one finds 

x - /^n fi(v) f / \ 3 v — arctan v 

e, { n, v ) = l +x , B{ n) + T J±l^, Mv) ^- — w —— j , (23) 

g,(n,v) = i + xjB ( n ) + nm+-v) ' /•(») = J, [-» + (! + » a ) . < 24 ) 

v _ _^_Q_ — cos ] 1 ^._ (25) 
c + 7 v ; 

These formulas are used for numerical calculations of the impedances. They have to be 
compared with the classical expressions in the local limit which follows from Eq. (fTTjl after 
the change to imaginary frequencies: 



z ° c = i z lT = v £( n) - ' £ W = 1+XlBin) + nWTVy (26) 



• Z i" — ^ -(12)= I + — ' 

The numerical result for Z s is shown in Fig. 01 as a function of Q for two values of Q. All 
calculations were performed for the parameters corresponding to gold at room temperature: 
7 = 3- 10~ 3 , vp = 1.4 • 10 8 cm/ 's, uj p = 1.37 • 10 16 rad/s. The impedance of the local theory 
is presented by the dashed lines. One can see that the nonlocal effect is very small for this 
polarization. The largest deviation from the local curves is just about of 2%. Obviously the 
s polarization cannot produce significant nonlocal correction to the Casimir force. 

A different situation is realized for p polarization, as shown in Fig. |S| One can see that 
there is a significant difference between the local and nonlocal cases. The deviation increases 
with frequency decrease and becomes larger for larger Q. This behavior has deep physical 
meaning, as explained below, and can appear only for the evanescent fields. Since in both 
cases the deviations from the local case are in the low frequency range, we analyze this limit 
analytically. 
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C. Low frequency behavior of impedances 



At low frequencies fi < 7, the variable v defined by Eq. (|25|) can be large if 7 < Vf/c « 
4.7 ■ 10~ 3 . Let us consider the impedances in the limit v ^> 1. In this limit the functions 
fi (v) and f t (v) in Eq. (J23I) and Eq. (J2U) behave as 

In the transverse dielectric function e t one can neglect 1 + %/b (fi) since the third term 
behaves as 1/fi at low frequencies. It gives for e t (fi, v) 

At: c 1 1 
3 Q cosh^fi 

For the the longitudinal function ej the phenomenological susceptibility xib (fi) again is 
negligible because it is responsible for the interband transitions at much higher frequencies 
but we cannot neglect the unit since the third term in (J23*|) does not depend on frequency 
at all and not necessarily large. For E\ (fi, v) one find 

E '< n -»> Bl+ <£d^) >=1+3 tes)'- < 29 > 

This expression describes the Thomas-Fermi screening of the longitudinal electric field. It 
has to be true [40] at fi < 7 and k much smaller than the Fermi wave number kp that is 
the applicability range of the Thomas- Fermi approximation. The latter condition k <^kp 
is also the condition for applicability of the Boltzmann equation. 

Substituting (J2SJ) and (j2U|) in Eqs. (j2Tj) and one finds for the impedances 

Z s = ~F{b), (30) 

y Q 1 1 ^ , 1 Q 2 v F Q 

Z P = 75 , — + -^GiP) « -7=-^— + ^G{b), 31 

n y/l + 3 (c/v F Q) 2 Q V3to c Q 

where the functions F (b) and G (6) are defined as 

00 2 00 • 2 

^ tc J ^cosh 3 x + b 3 ' ^ 7r J "^cosh 3 x + 6 3 ' ^ ^ 
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Q V 4 v F J 

The functions F(b) and G(b) can be found explicitly but the result is cumbersome and 
inconvenient for analysis. For this reason we calculated the integrals in Eq. (|32|) numerically 
presenting explicitly only the asymptotics at b <C 1 and b ^> 1. The functions F(b) and G(b) 
are shown in Fig. El The asymptotic behavior of these functions is 



F(b) 



±b\ b«l 



4 1 

3^ b 



+ ^3 (In 2b -1/2), 6>1 







G(6) = { 2 157r ' (34) 

aTil-i ( ln2fe + 1/2), 6>1 
The known result for the Leontovich impedance for the strong anomalous skin effect 
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42} is easily reproduced if we take in the equations above the limit Q — > 0. In 
this limit the parameter b goes to infinity and the contribution of the transverse dielectric 
function is the same for both polarization: F(b) = G(b) = 4/3\/36- The contribution from 
ei in Z p disappears in the limit Q — > 0. Hence, the impedances will coincide with each other 
and they are given by the classical expression for the strong anomalous skin effect continued 
to imaginary frequencies 

z s (o, n) = z p (0, n) = z (n) = ^ ^ fi2 ) " ■ (35) 

However, if Q is nonzero there is a small enough frequency where b is not large anymore and 
Eq. ()35|) is not applicable. When Q is so small that b 1 the impedance Z s approaches 
the limit Q/Q. The same limit is realized for the local impedance Z l ° c in Eq. (J26|) at very 
low frequencies when one can neglect e (Q) in comparison with Q 2 /Q 2 . 

For p polarization at nonzero Q the contribution of St in the impedance decreases with 
Q but the contribution of E\ increases as I/O, (see Eq. ()3ip ) and dominates in Z p at low 
frequencies. It is in agreement with our numerical calculations in Fig. |S| Indeed this is the 
result of the Thomas-Fermi screening. The same effect is not realized for the propagating 
fields. In this case the ratio Q/Q = sini? < 1 is restricted. Since vp/c = 4.7 ■ 10~ 3 is small, 
the variable u (|16|) is small nearly everywhere in the integration range and the function 
fi (u) ~ 1. Therefore the local limit is realized instead for the longitudinal contribution in 



Z p . 
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The behavior of the impedance for the s polarization at low frequencies is a sensitive 
matter for the temperature correction to the Casimir force. One of us (VBS) in collaboration 
with M. Lokhanin analyzed this problem with the Leontovich impedance (J35|) . As follows 
from the discussion above this analysis has to be reconsidered taking into account not only 
different behavior of Z s at very low frequencies but also significant deviation of Z v from the 
local impedance in this range. 



III. NONLOCAL CORRECTION TO THE CASIMIR FORCE 

In this section we are going to estimate the correction to the Casimir force due to non- 
local effects at frequencies smaller than uj p . The restriction on frequency is connected with 
the use of the Boltzmann approximation for the dielectric functions (|2*3"|) . (J2IJ). In this ap- 
proximation we cannot describe the plasmon excitations. Of course, one could use more 
general dielectric functions like those in the self-consistent-field approximation j^j] to an- 
alyze all the nonlocal effects. However, we think it is reasonable to separate the effects of 
different physical origin. Influence of the plasmon excitations on the Casimir force has been 
already evaluated |38J] using the hydrodynamic approximation for the longitudinal dielectric 
function, but the correction to the force due to the anomalous skin effect never has been 
calculated before. Only specific questions concerning the temperature correction have been 
addressed in the literature 3]. By anomalous skin effect we refer not only to the strong 
anomalous skin effect that is realized when the electron mean free path is larger than the 
field penetration depth, but to all the nonlocal effects that happen at frequencies smaller 
than Up. 

We will consider only the force in the zero temperature limit. Thus, the Casimir force will 
be calculated without the temperature correction but all the other parameters characterizing 
the material, especially the relaxation frequency u T , will be kept at finite temperature. 
The results of the previous section for the impedances Z s and Z p are important for the 
temperature correction problem but this question will be considered elsewhere. 

It is known that when a metal is described by the surface impedances, the Lifshitz formula 
for the Casimir force remains essentially the same [25, 38] as when the metal is described 
by the local dielectric function. Only the reflection coefficients have to be expressed via the 
impedances. At nonzero temperature the Lifshitz formula includes summation over the 
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Matsubara frequencies ( n , defined for our dimensionless frequency as 



_ Cn_ 2%nkT 

" "77. 



(36) 

To get the Casimir force at T = we have to integrate over the continuous variable Q. In 
the dimensionless variables Q and Q the Casimir force between two metallic plates separated 
by the distance a at T = is 



F pp (a) 



he 



oo oo 



27T 2 5 4 



dn / dQQ^Q? + Q 2 



rj 2 exp (2cVfi 2 + Q 2 ) - l) + 







where 



r p 2 exp (2dy/W + Q 2 ) - 1 



(37) 



a c 
d — —, 5 = — « 21.88 nm. 



5' 



(38) 



Here 5 is the penetration depth for gold. The reflection coefficients follows from Eq. (JTSJl 
after continuation to imaginary frequencies 



n - Vfi 2 + Q 2 Z S (fl, Q) 



'" Q + Vfi 2 + Q 2 Z S (Q, Q) ' ^ y/fi 2 + Q 2 + (O, Q) 
The impedances Z S)P are calculated according to Eqs. (|2*T|) - (|2*3|) . 

The force between a sphere and a plane can be calculated with the help of the proximity 
force approximation 



^n 2 + Q 2 - VZp (O, Q) 



(39) 



^ (a) 



5lj which gives the following expression 
hcR 



oo oo 



J dttj dQQ In (l -r s 2 exp (-2dy/n 2 + Q 2 )) + 



In ( 1 - r; 



;ex P (-2dy/n 2 + g 2 



(40) 



where i? is the radius of the sphere. 



A. Numerical procedure 



First we calculate the force between parallel plates F pp e (a) in the local limit with the 
Drude dielectric function and local impedances (j^Sj) - The actual calculations were performed 
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'or the dimensionless relaxation frequency (see ([20))) 7 = 4- 1CT 3 . This value is the best fit 
17 1 of the handbook optical data for gold at low frequencies. In Fig. [7|we show the 

reduction factor 77(a) defined as the ratio of the calculated force to the original Casimir force 

(HJ); this is 

Via) = F f^-. (41) 
F c {a) 

The force calculated using the Drude model (dashed line) is smaller than that calculated 
using the handbook optical data for gold (solid line). The solid line coincides with the 
reduction factor given in Ref. 

The nonlocal correction is calculated without the empirical susceptibility xib (<^) intro- 
duced in Eqs. (|12|) and (|13|) so we have to remember that the relative nonlocal correction 
will be smaller than the calculated one on the value of the order of the relative difference 
between the curves in Fig. [7| (8% at a = 100 nm). 

The force calculation in the nonlocal case is quite complicated because one has to make 
three integrations with high precision, one to calculate the impedances and two to calculate 
the force. It is much more easy to calculate not the force itself but integrate the difference 
between nonlocal and local integrands. In this case there is no need to perform high precision 
calculation of the integrals since we have to know the correction due to nonlocality with the 
precision of about of 10%. Actual calculation of the difference 

SF„ (a) = (a) - F l ;;{a) (42) 

were made with the relative accuracy of 1%, while the impedances ()21|) . (|22~j) were calculated 
with the relative precision of 10~ 6 . The integrands for s and p polarizations defined as 



Sf s , P = Qy/W + Q* [(rjexp {2dy/fF+&) - l) * - (r., - r l ° c p ) 



(43) 



are presented for a = 275 nm in Fig. Ufa) and Fig. (Hfb), respectively. Both of them are 
negative as it should be, since the force decreases due to the nonlocal effects. It is interesting 
to notice that Sf s is nonzero in a very narrow range of small Q. In contrast, the integrand 
for p polarization 5f p is nonzero in a broader range of Q (pay attention on different scales in 
figures (a) and (b)). Nonlocal effects are always significant in a wider range of Q < 1. With 
the decrease of separation a the integrand for p polarization increases in the magnitude and 
becomes wider in both directions Q and Q. The integrand for s polarization decreases in 
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magnitude and widens only in Q direction. Thus, the contribution of s polarization in the 
force correction is always smaller than that for the p polarization. 

The results for the relative correction 5F PP (a) / F pp rude (a) due to the nonlocal effects are 
presented in Fig. El The solid line gives the resulting correction, while the dashed and 
the dotted lines represent the contributions of p and s polarizations, respectively. One can 
see that the correction is small but not negligible. Contribution of s polarization increases 
when 7 becomes smaller, but even for 7 = 10~ 5 this contribution is still on the level of 
0.2%. We can see that large deviation of the impedance for p polarization from the local one 
that happens at low frequencies is not very significant for the force. This is because in the 
reflection coefficient r p the impedance enter as QZ p so that 1/fl behavior of Z p is suppressed 
in the reflection coefficient. 

Similar calculations were made for the sphere-plate geometry. The relative correction 
SF sp (a) I F® rude (a) together with the separate contributions of p and s polarizations is 
shown in Fig. The behavior of the curves is quite similar to that for the plate-plate 
geometry. Only the absolute magnitude of the correction is smaller. 



IV. DISCUSSION 

The theory described in Sec. |H] provides a solid ground for the impedance approach in 
the Casimir force calculation. Specifically it allows correctly to take into account nonlocal 
connection between the displacement and electric fields. In this paper we restricted ourselves 
by the nonlocal effects happening at frequencies smaller than oo p . This restriction is due to 
used Boltzmann approximation for the nonlocal dielectric functions (fT2"|) - (jl6|) . However 
the equations for the impedances Q, (jlOJ) are much more general. For specular electron 
reflection off the surface these equations are true for arbitrary dielectric functions e% (k, 00) 
and ei(h,u) with the only condition that these functions exist. Therefore, all nonlocal 
effects can be described on the same basis. In particular, for metals the most general 
dielectric functions for a free-electron gas were found in the self-consistent-field (or Lindhard) 
approximation with the necessary modifications to include a finite relaxation time 40]. In 
this approximation et(k,uo) has the following form 

S t (fc, U)=l — rf t [U, Z) , 

00 [00 + %U0 T ) 
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where 



f t (u,z) = -(z 2 + 3u- 2 + l) 



32z 



1"( 



-1\2 

: — u ) 



2 , fz-U^ + l 

m 



z — u 



-i 



i\ 2 



1 - (z + u- 1 ) 



111 



z + u 



(44) 



+ - 1 

Here u defined as before by Eq. (JTfij) and z is z = k/2k F . The longitudinal dielectric function 
has a little bit more complicated form: 

ei{k,u) = l + [e w -l) 1+t In 

io 2u oo VI — u 



£w = 1 + k^ fl ( u > ^ 



In 



z - u~ l + 1 



z — u 



-1 



l-^+^i- 1 ) 2 



111 



z + u' 1 + 1 



(45) 



z + vr x - 1 

All the other approximations for the free-electron gas can be found from f)44j) . (|43|) in definite 
limit cases. For example, the Boltzmann approximation (jl2|) - (fTT)j) follows from ([44)1 . ()45|) 
in the limit z — > 0. These expressions for the dielectric functions allow to perform detailed 
investigation of the high frequency region oj ~> u p which gives more significant contribution 
in the Casimir force due to excitation of the propagating charge density waves in the metal 



We considered here only specular electron reflection off the metal surface. It is justified for 
the AFM experiment where the root mean square (rms) roughness of the surface (1 nm) 
was much smaller than the mean free path (30 nm). However, in the MEMS experiments 

0, the rms roughness was comparable with the mean free path and approximation 
of specular reflection fails. In this case the diffuse reflection of electrons off the surface is 
more suitable. For the diffuse reflection the impedances are not represented by the Eqs . 
(0), anymore. Instead one has to use the impedances for the diffuse reflection 4l|. 
There is no problem with Z s which is expressed via e t (k, oj) but situation with Z p is much 
more complicated. This occurs because of the destruction of translational invariance in 
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the direction normal to the surface 



4l| . Although it is possible to calculate both of the 



impedances in the diffuse case, we do not think it is reasonable to do for the anomalous skin 
effect. This is because the nonlocal correction is smaller than the uncertainty in the Casimir 
force due to the roughness. The roughness correction to the force is usually evaluated using 
the proximity force approximation (see, for example, [9|). Recently it has been pointed out 
52I ] that this approach is valid only for long wavelength deformations of the plates. The real 
surfaces of deposited gold films always roughed on very different scales |53| and the short 
wavelengths will bring uncertainty in the estimate of the force. 

The impedances (J0J), (jlOj) with the nonlocal dielectric functions (|12j) - (|lfjjl are well known 
in the optics of metals but in this paper we considered them in the near field range where 
the nonlocal effects were unexplored. In this sense the Casimir force is a unique problem. 
Significant contribution in the force comes from fluctuating fields in the near field region. 
Since the force has to be predicted with high precision, it is important to take into account the 
nonlocal effects. Though we have found here that the anomalous skin effect gave observable 
but minor correction to the force, the other nonlocal effects, such as plasmon excitation, can 
give more significant correction. This paper just provide a regular basis for the calculations 
of this kind. 



V. CONCLUSIONS 



A complete calculation of the Casimir force that can be accurately compared with exper- 
iments requires a full optical characterization of the involved materials. This is complicated 
due to the various factors that can modify the optical properties. In this work we described a 
systematic way to take into account the nonlocal effects in the material. It was stressed that, 
in general, a metal had to be described with two different surface impedances corresponding 
to s and p polarizations and these impedances depend not only on frequency but also on 
the wave vector along the metal surface. As a specific problem we considered the correction 
to the Casimir force in the region of the anomalous skin effect {u < uo p ). This region is 
characterized by the nonlocal dielectric functions (longitudinal and transverse) that can be 
obtained in the Boltzman approximation. The impedances are completely defined by these 
functions. 

It was demonstrated that the exact impedances were differed from the approximate Leon- 



20 



tovich impedance. The latter one caused confusion in the literature so our analysis resolved 
the problem and gave a proper description of the impedance approach in the Casimir force 
calculation. It was emphasized that the significant contribution in the force came from the 
evanescent fields. For these fields the impedances can be found by the analytic continuation 
and the procedure is well defined. The contribution of the nonlocal effects in the impedances 
was found to be quite different for propagating and evanescent fields. Specifically for the 
evanescent fields the impedance for p polarization deviates significantly from the local one 
that is the result of the Thomas-Fermi screening. For s polarization the nonlocal contribu- 
tion in the impedance is more significant for the propagating fields than for the evanescent 
ones. 

In the impedance approach the Casimir force can be found from the same Lifshitz formula 
in which the reflection coefficients are expressed via the impedances. We calculated the 
nonlocal correction to the force in the region of anomalous skin effect at zero temperature. In 
spite of significant deviation of Z p from the local impedance the nonlocal reflection coefficient 
r p deviates from the local one only slightly. For the s polarization the effect is even smaller. 
For this reason the total contribution of the nonlocal effects in the Casimir force is on the 
level of 0.5% at small separations. It is a minor effect within the levels of detectability of 
present experiments, but smaller than the corrections introduced due to sample roughness. 

We did not considered in this paper the temperature correction though it is clear from 
the analysis of impedances that anomalous skin effect will be important for the temperature 
correction. A new phenomenon observed here is that the reflection coefficient r p for p 
polarization is not going to 1 in the zero frequency limit u — > 0. This behavior is the result 
of the Thomas-Fermi screening. 

The technic developed in this paper can be applied to calculate the contribution of the 
other nonlocal effects such as plasmon excitation at uj ~ u p . These effects are expected to 
give more significant correction to the Casimir force. 
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FIG. 1: Choice of the coordinate system for the incoming wave. The angle of incidence is given by 
i?. The electric fields in the p polarized wave E p and s polarized wave E s are shown. 
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FIG. 2: Absorptance as a function of the dimensionless frequency Vt at normal incidence = 0. 
Local and nonlocal cases are represented by the dashed and solid lines, respectively. At $ = there 



is no difference between s and p polarizations. Parameters were chosen as in Ref. 
v F = 0.85 • 10 8 cm/s. 
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FIG. 3: Absorptance as a function of Q for p polarization at the incidence angle 1? = 75°. Local 
and nonlocal cases are represented by the dashed and solid lines, respectively. Parameters were 
chosen as in Fig. |2J 
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FIG. 4: Numerically calculated impedance Z s as a function of dimensionless frequency f2 for 
two values of the dimensionless wave numbers along the plate Q. Solid line describes nonlocal 
calculations; the dashed line present the local case. Maximal deviation between local and nonlocal 
curves is about 2%. Gold parameters were used for calculation: 7 = 3- 10~ 3 , vf = 1-4 • 10 s cm/ s. 
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FIG. 5: Numerically calculated impedance Z p as a function of frequency for two values of the 
wave numbers along the plate Q. Nonlocal and local cases are shown by the solid and dashed lines, 
respectively. 




FIG. 6: Numerically calculated functions F(b) and G(b). 
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FIG. 7: The reduction factor F pp (a)/F c (a) in the local case as a function of the separation a 
calculated with the handbook data (solid line) and with the Drude model for the dielectric function 
(dashed line). 




FIG. 8: Integrands for s polarization (a) and for p polarization (b). Note different scales in Q axes. 
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FIG. 9: The relative correction to the force due to nonlocal effects for the plate-plate geometry. 
The solid line presents the resulting correction. The dashed line gives the contribution of the p 
polarization and the dotted line gives the contribution of the s polarization. Gold parameters were 
used for calculations. 
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FIG. 10: Same as Fig. 03 but for the sphere-plate geometry. 
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